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Abstract 

We study the distribution function, P{p), of the effective resistance, p, in two and three- 
dimensional random resistor network of hnear size L in the hopping percolation model. In this 
model each bond has a conductivity taken from an exponential form a oc exp(— «;r), where k is a 
measure of disorder, and r is a random number, < r < 1. We find that in both the usual strong 
disorder regime L/k'^ > 1 (not sensitive to removal of any single bond) and the extreme disorder 
regime L/k"^ < 1 (very sensitive to such a removal) the distribution depends only on L/k'^ and 
can be well approximated by a log-normal function with dispersion bn'^ /L, where 6 is a coefficient 
which depends on the type of the lattice. 
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INTRODUCTION 



The concepts and methods of percolation theory are widely used to explain many phe- 
nomena in physics, classical as well as quantum (for review see e.g. Refs. 

1133). The 

canonical model for studying transport properties of disordered systems is percolation on a 
lattice. Usually it is also assumed that the conductivity between neighboring lattice sites 
may be defined as either finite, or zero (i.e., either conducting or insulating) without loss 
of generality. This model [which we denote as the bond (or site) percolation] has been ex- 
tensively studied and is understood quite well. For the description of the nearest neighbor 
hopping in granular materials it is much more natural to define the conductivity between 
two neighboring lattice sites (labeled as "i" and "j") by 0"^ oc exp[— rj^/ro — eij/ksT], where 
Vij is the distance between the two sites, vq is the scale over which the wave-function outside 
the grain decays, e^j is the energy difference between grains, and T is the temperature. Here 
we neglect the thermal hopping term (high temperature regime) and consider only nearest 
neighbor hopping. This behavior may be captured by a lattice model for which jj, ^,0,0,0 

aij = aoexp[~Kr{ij)], (1) 

where k is a measure of disorder, r{ij) is a random number taken from uniform distribution 
in the range (0,1), and do is a dimension coefficient P, |^. We shall name this model the 
hopping percolation model. 

One might expect that such small differences in the formulation of the disorder in the 
aij (namely, (Xij = 0, 1 or aij = (Tq exp[— Kr(ij)]) will lead to no important difference in the 
global conductance properties of these systems. Quite surprisingly, recent experiments on 
the conductance of granular material might suggest otherwise. Specifically, the number of 
red bonds (which are critical for current) expected in framework of the traditional percolation 
theory is proportional to L^/'^, where u is the percolation correlation critical exponent and L 
is the system size [l^. Thus, for a typical experimental set-up of 10^ grains, one expects 10^ 
red bonds. On the other hand, measurements of transport through Ni granular ferromagnets, 
indicates a much lower number of red bonds (typically of order one) 0]. In a recent paper 

we have attributed this difference to the fact that the estimation of the number of red 
bonds (L^^'^) is based on the bond percolation model, while for the hopping percolation we 
expect a transition to a regime of extreme strong disorder in which a single red bond governs 
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a. b. 

FIG. 1: (a) A square bond percolation net of resistors with random resistivity given by Eq. 
where (ij) denotes the bond between sites i and j. (b) A site percolating network. The resistivity 
of all four resistors within a dashed circle is determined by a single random number r(ij), where 
ij denotes the labeling of the grid point 

the behavior of the system. The onset of this regime scales as to k/L^^'^. 

It is important to note that in contrast to the traditional bond (or site) percolation 
model, in which the system is either a metal or an insulator, for the hopping percolation 
model the system always conducts some current. Hopping conductivity (i.e., exponential 
local resistance) is always associated with strong disorder. As was shown in Ref. , there 
are two regimes within this strong disorder: a regime which is not sensitive to the removal 
of a single bond, as expected from the usual percolation theory, termed the usual strong 
disorder regime 11 . While for even stronger disorder a regime which is very sensitive to the 
removal of a specific single bond exists, denote as extreme disorder. In the extreme disorder 
regime a single bond can determine the transport properties of the entire macroscopic system 

m 

The remainder of this paper is arranged as follows. In Section „ we describe our model 
and the numerical approach. In Section „ we present some numerical results, followed by a 
brief discussion in Section „ 

MODEL 

We perform large-scale Monte Carlo simulations for calculating transport in these sys- 
tems. We build a bond-percolating Miller- Abrahams like resistor network (see Fig. Q and 



3 




FIG. 2: A semi-log plot of the averaged resistance pe vs. k. In the case of site percolation 
{pc = 0.592746) the slope of the curve is close to 0.6, while for the case of bond percolation 
{pc = 0.5) this slope is equal to 0.5 [cf. with Eq. ((21)]. The system sizes shown are: L = 20 (solid 
circles), 40 (open squeres), 100 (open stars). 
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IJ), but assume the conductivity of each resistor to have the form given in 
Eq. (P). Then solve the corresponding set of hnear Kirchhoff equations and calculate the 
total effective resistance pe for two dimensional (2D) and three dimensional (3D) networks 



We b^in by calculating the average effective conductivity (Xe. The approximate expres- 
sion P, UQ| for the effective conductivity cTg, of a random resistor network with local 
conductivities given by Eq. in 2D is 

ae = a,e-''^\ (2) 

In Fig. 121 we show this dependence (in terms of resistivity pe = l/cTg) for both site and bond 
percolations (see Fig. for different lattice sizes. In Ref. jsl it was shown that in the 
limit K ^ oo, Eq. ^ is exact. It is easy to show that in the case of 2D random resistor 
bond network (for which pc = 0.5), Eg. ((21) follows immediately from the Keller-Dykhne 
theorem and is exact for arbitrary k [l^. Similarly this results can be found in framework of 
the symmetric self-consistency effective-medium approximation (EMA) rewritten for many 
component composite [l^. From Eq. ((21) follows that the effective conductivity de depends 
on K and does not depend on the system size L. For finite L, Eq. ((2} represents the mean 



leDercolation correlation lenffth ^ oc (p~Pc) " 

"inn rf ' 

Ij, y, 1^ ) • It was also shown 5} that h describes 



conductivity over all configurations of the disordered system. 

Next we study the fiuctuations of the resistance, p = l/a, from the mean value pe = V^e 
for individual systems of finite size L. We perform numerical calculations of the probability 
distribution function P{p) (i.e., the probability that the total resistance of the system is p) 
as well as the variance var(p) as a function of L and k. As shown in Ref. |5j, the relative 
variance (in contrast to pe), strongly depends on L and k only through the scaled variable 
h = L/ . Here v is the critical exponent of t 
(in 2D z/ = 4/3 ~ 1.33, while in 3D z/ ~ 0.88 
the transition from strong disorder [h > 1) to extreme disorder {h < 1). 

RESULTS 

Here we present numerical results suggesting that P{p) also depends only on L/ . In 
order to verify and to quantify this hypothesis, we study numerically P{p) for systems of 
different sizes L and different disorder k, but with the same value of h (see also Ref. jl3|)- 
In Figs. El and m we show P{p) vs. p for the cases of strong and extreme disorder. All data 
corresponding to the same parameter h scale according to the same law [see Figs. Efc) and 
Efg)]. Thus, our results suggest that P{p) is a function of both p/pe and /i, i.e., 

P{p) = -f{^-A- (3) 



Pe \Pe 

Here h determines the form of the function and -P(p) for a fixed h depends only on p/pe- 
Fig. Ol^c) suggests that in the strong disorder, h controls the width or standard deviation 
of the rescaled distribution. Since the standard deviation increases when h decreases, we 
assume that the standard deviation 6 ~ bpe/h = bpen'^/L, where 6 is a parameter which 
depends on the type of lattice. Indeed, when we plot in Fig. Efd) P{p/pe)S vs. (p — Pe)/5, 
a collapse of the two plots shown in Fig. ^c) is obtained. The functional form obtained in 
Fig. Efd) suggests that the probability distribution can be approximated by a Gaussian 

P(p) ~ iV2^6)-' exp [- (p - Pe)' /26'] . (4) 

Indeed, the dashed line in Fig. Efd) represents a good fit to the Gaussian given by Eq. 
However, Eq. (0} can not approximate the asymmetric form of -P(p) at extreme disorder [see 
Figs. Efg) andO^h)]. We suggest, as will be justified below, that -P(p) can be approximated 
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FIG. 3: (a)-(d) The probability distribution for the case of usual strong disorder (h > 1). (a) A 
typical form of P{p) vs. p for h = 8.1. (b) Similar to (a), but shown in log-log scale for L = 300, 
K= 15, h = 8.1 (o), L = 500, K = 22, h = 8.1 (a), L = 100, k = 10, h = 4.7 (*), and L = 200, 
K = 16.84, h = 4.7 (•). (c) Scaling plots peP{p/pe) vs. p/pe for the results shown in (b). Only the 
data with the same value of h scale to a unique function, (d) Scaling plots of peP{p/ Pe)/{L/ k^) 
vs. {L/ k'^){p/ P(. — 1) lead to the collapse of the data shown in (b) and (c). The dashed line in 
(d) is the Gaussian distribution Q with h = 0.2. (e)-(h) Similar to (a)-(d), but for systems with 
extreme disorder {h < 1): L = 50, k = 30, h = 0.54 (A), L = 30, k = 20 A, h = 0.54 (■), L = 40, 
K = 40, h = 0.3 (□) and L = 20, k = 23.78, h = 0.3 (x). (h) Similar process like (d), does not 
lead to collapse in the case of extreme disorder. The lines connecting the points for h = 4.7 in (c) 
and for h = 0.3 in (g) and (h) are guides to the eye. 



(in all regimes of disorder) by the log-normal form 



Pip) 



exp 



ln^(p/Pe) 



(5) 



where fi = 6 / pe = hn^ /L. In fact, Eq. © includes also the usual strong disorder case, since 
in the latter case ln^(p/pe) — (p/Pe — 1)^ and Eq. ^ reduces to the Gaussian form (jlj), 
while at extreme disorder (/i ^ 1) the exponent function in Eq. (0) tends to 1, and P{p) 
transforms to the power-like dependence ~ 1/p (see Fig. E}. 
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FIG. 4: A log-log plot of P(^p) vs. p. By increasing disorder the log-normal distribution transforms 
to the power-law (a) L = 50, k = 30 (/i = 0.54). (b) L = 20, k = 35 (/i = 0.18). Inset: The 
exponent s (of the power law p^**) vs. h. By increasing disorder (i.e., decreasing K) the exponent 
s tends to 1. L = 50, k = 30 (/i = 0.54); L = 14, k = 18.19 (/i = 0.3); L = 20, k = 35 (/i = 0.18); 
L = 14, «; = 35 (/i = 0.12). v = 4/3. 

From Eq. (0) it follows that at /i ^ 0, the distribution function P(p) reduces to a delta- 



function: limP(p) 



lim 



ln(p/pe)^ 



I.e., K 



P '-'o V2n,, - l^^^"^ p -In p^) = 6{p-pe). Therefore, at p ^ 

or L ^ oo) the total resistivity of the system is exactly poe^"'^ and has no size 
dependence: lim J pP{p)dp = p^. 

It should be noted that a log-normal distribution of resistances is found in quantum 
models of hopping conductivities (see e.g., Ref. and references therein), while here it 
is demonstrated for classical exponential disorder ([T]). Moreover, our result (0) yield the 
specific analytical form of P{p), which includes the dependence on k and L for all regimes 
of disorder. 

In Fig. Efa) we test Eq. (0) by comparing it to simulation results. It is shown that the 
numerical results of the 2D resistance pP{p/ pe) / {L/ k") scale vs. (p/pe)^^'^", as predicted 
by Eq. (0) for both strong and extreme disorder. A similar plot is presented in Fig. Efb) for 
a 3D lattice. Although for the 3D case Eq. Q is not exact (since Keller-Dykhne theorem 
exists only in 2D), nevertheless the approximated expression cTg ~ aQu'^e'^"'^ is known 0, 
0], resulting in the distribution law Since in Eq. Q the parameter n appears with 
the prefactor pc, we should expect that pc enters into the parameter p of Eq. © as /i = 
a{pci^Y / L. Comparing the values h = 0.2 observed for the square bond percolation lattice 
{pc = 0.5) and 6 = 0.18 for the cubic site percolation {pc = 0.31 
The dependence of p on pc is in agreement with result of Ref. 



6), we find that a = 0.503. 
l^ . These results strongly 
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FIG. 5: (a) A scaling plot of pP{p/pe) vs. {pjp^^l'^'' of the data plotted in Fig. 01 (2D case). 
The dashed line represents the analytical result @, with h = 0.2 and = 4/3 for 9 systems with 
L = 500, K = 22 (A); L = 300, k = 15 (■); L = 200, k = 10 (•); L = 100, k = W (*); L = 50, 
K = 30 (A); L = 40, K = 40 (□); L = 30, k = 20 (o); L = 20, k = 40 (x); L = 20, k = 30 (0); 
(b) Similar to (a), but for 3D, with b = 0.18 and ly = 0.88 for three systems with L = 20, k = 15 
(□); L = 26, K = 15 (A); L = 10, k = 6.8 (*). (h = 1.85). The dashed line is the analytical result, 
Eq. ©. (c) ^{pc) = ^[ln(/5)/'t] vs. Pc, where Pc = lii(p) for different values of the ratios L/k'^: 
L = 100, K = 10, /i = 4.68 (*); L = 60, k = 30, /i = 0.65 (o); L = 40, k = 20, /i = 0.74 (□); 
L = 20, K = 30, /i = 0.22 (A); (d) Scaling ^\n{p) / n]/ {L / kJ"-^) vs. [ln(p)/K - 0.5](L/k'^"1) with 
V = 4/3. The values of the ratios L/k^ are the same as in (c). The dashed line is the analytical 
fit, as derived from Eq. ©. 

support our proposition that Eq. © describes well the distribution in all ranges of disorder. 

The variance var(p) can be expressed as (p^) — Pg, where (p") = p^P{p)dp, and 

Pe = (p). For large enough k, (p") = p'^P{p)dp = p^e^^'^P'^+i^^i^^ and the relative 
variance takes the form 

[var(p)]^/VPe = [e'^'(e'^'-l)]'/' (6) 
(see Ref. j3|). Fig. IHl presents numerical results showing that the relative variance scales as 
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FIG. 6: (a) A semi-log plot of the relative variance [var(p)]-'^/^/pe vs. k for various sizes of the 
system: L = 10, 14, 20, 30, 40, 50, 100, 200, 300 (from top to bottom), (b) A semi-log scaling plot of 
the same quantity vs. = bK^ /L, where b = 0.2, v = 4/3. The dashed line represents Eq. 0. 

a function of /i = bn'^/L in accordance with Eq. 

Next we shall present analytical arguments for the log-normal distribution (0). According 
to the central limit theorem |2(J], if the values of In p are normally distributed, then the values 
of p should follow the log-normal distribution. Assuming In p = upc [see Eq. (j21)] for all p, the 
distribution P{k,~^ In p) is simply the distribution of the percolation threshold $(pc) which 



is normally distributed (e.g., Refs. 2l|). Indeed, in Fig. El^c) we show that <l'[K~^ln(p)] 
approximately follows a normal distribution centered at pc = 0.5. Thus, the distribution of 
p should be log-normal as in Eq. (jSJ. 

Using the above assumption Inp = p/n = upc, it is possible to evaluate the distribution 
$(pc) and its standard deviation 6p^. One can write a simple relation ^{pc)dpc = ^'{y)dy, 
where y = npc, and get $(pc) = ^'{u)'^ = f^^iv)- Therefore, $(pc) = K$'(lnp) ~ 
exp[—K^{pc — Pc)^/2/i^] = 6p^ exp[—{pc — Pc)^/25pJ, where Pc is the mean value of 
the percolation threshold and 6p^ = p/n is the standard deviation of $(pc)- This form of 
$(pc) is supported by our numerical simulations shown in Figs. |31|^c) andEfd). 

This specific form for 6p^ = bn'^^^f L in the hopping percolation model should be compared 



to 6p^ = L known for the bond percolation model 10l|21|. This furt 



differences between the bond percolation model considered in Ref. 
percolation model considered here. 



ler emphasizes the 



2l| and the hopping 



SUMMARY 



In summary, we find the specific form of the resistance distribution in the hopping per- 
colation model. For all ranges of strong disorder /t and lattice sizes L, the distribution is 
log-normal and depends only on the ratio k'^/L, where v is the correlation exponent for the 
bond percolation case. Assuming the relation p = exp(/tpc) for finite systems leads to a vari- 
ance oi Pr, Sp^ = /L, which is different from 5p^ = L^^^^ known for the bond percolation 
model j2l|. Our results may be relevant to ac conductivity measurements in such systems. 
By appropriate choice of frequency one can detect regions of size smaller than k'^, where a 
crossover in behavior from extreme to usual strong disorder behavior is expected. 
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Ministry of Absorption of the State of Israel. 
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